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During the first order phase transitions the objects of a new phase con- 
sume a mother phase which leads to some density gaps near existing embryos 
of a new phase. This fact is rather obvious but complete theoretical de- 
scription was given only in 0. Now this fact in widely recognized but still 
there are some attempts Q to reconsider this approach which is now known 
as the model of depletion zones ||. Here we shall analyze || and show the 
useless of reconsideration proposed in J|]. 

We don't consider here the picture of the gobal evolution of a system - in 
M it is so primitive that can not give any reliable quantitative results but only 
some very rough estimates without any justifications. We needn't to do it 
because this has been already done in |2| with all necessary justifications. So, 
we consider only the profile around the solitary embryo, which was the matter 
of interest in M. Why it is possible to consider only the solitary embryo? 
An answer has been given in |J, in Q this question remains unsolved. Due 
to results of we shall consider namely the solitary embryo. 

We assume that the substance exchange regime is a diffusion one. Then 
there is a gap of a vapor density n near an embryo. Then the density n differs 
from a density n(oo). This difference initiates a variation in a nucleation rate 
/ and causes the effect of depletion zones. So, it is very important ro get for 
n(r) — n(oo) (r is a distance from a center of an embryo) a correct expression. 
This expression is a crucial point in consideration of an effect of depletion 
zones in nucleation kinetics. Certainly, it has to be a function of time t. 

In |4|] the following expression is proposed 

n(r,t) = n s (r) + (n(oo)+n(0))%(^£) (1) 
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/ \ /x n(oo) — n(0) „ 

njr) = n(oo) — — i? 2 

r 

Here n(0) is equilibrium concentration on the surface of nucleus, D is a 
diffusion coefficient, R is the nucleus radius, is the Laplace probability 
integral. 

Here the sign + was used instead of + in the second term of ([!]). Correct 
solution in the neigbourhood of an embryo is || 

n — n(oo) R r — R , . 2 



— er f c C -7= ), er f c ^ = ~f^l z ex p( - £ 2 )*; ( 3 ) 



n(0)-n(oo) r 2y/Dt " 

But even with this solution one can not construct the correct form of the 
nucleation rate profile. We shall outline this point. 

There are following necessary suppositions to get this exact solution for 
concentration profile: 

• The radius of an embryo is constant in time 

• The boundary condition at r = R is constant in time. 

Both assumptions are rather approximate. As a result the above solution 
is valid only near the surface of an embryo. It is easy to see from the following 
arguments. Let 5t be the characteristic time when the relative variation of 
intensity of vapor consumption v is already essential 

\v(t + 6t) -v(t)\ ~ v(t) 

Then it follows that only the distances 

r - R ~ AD5t = 5R 

or smaller ones can be considered on the base of solution (||). At all distances 
r — R > 5R we have to take into account that earlier the intensity of vapor 
consumption was smaller. Really 

v ~ du/dt 

where v - a number of molecules inside the embryo. As it is known 

du/dt = const z/ 1/3 ~ t 1/2 
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and we see that the intensity isn't constant in time. 

It seems that the small distances are the main ones. Certainly the gap is 
greater the smaller the distance is. But the small distances aren't essential 
because the functional form of the nucleation rate can be approximated 

as 

I if) = I(r = oo) exp(r(n(r) — ra(oo))/n(oo)) (4) 

where T is a big parameter. It is approximately equal to v. Note, that the 
sence of nucleation rate as a probability to appear for an embryo was used 
here (details see in 0). So, the distances where 

—T(n(r) — n(oo))/n(oo) ~ 1 
are the most interesting. When 

\T(n(r) — n(oo))/n(oo)| <C 1 
we have / = /(oo) and the is no gap of concentration. When 

— r(n(r) — n(oo)) / 'n(oo) 1 

we have 1 = and the is no nucleation. So, this region isn't interesting. 

The distances T(n(r) — n(oo))/n(oo) ~ 1 because of T ~ v — ► oo corre- 
sponds to very small n(r) — n(oo) and, thus, to big values of distances ~ 5r 
from the embryo. Ordinary 5r ^> 5R and there is no possibility to use (^). 

Instead of (|3]) in [|] another approach was used. This approach is based 
on the strong inequality 

5r > R (5) 

Why this inequality takes place? When we consider classical nucleation, 
i.e. transition of a supersaturated vapor into a liquid phase everything is 
clear. We have a strong inequality 

v v /v t > 1 

where v v is a partial volume for one molecule in a vapor, vi is a partial 
molecule in a embryo. This inequality makes obvious that the final vol- 
ume after the whole process of condensation is small. As an example of 
condensation of water vapor in normal external conditions. Suppose that 
supersaturation is somewhere ~ 5. Then because v v /v[ ~ 1000 we see that 
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the final volume of a new phase is very small ~ 1/200 and the characteristic 
distance 5r ~ (200) V 3 R > R. 

We see that the value 5r has to take the same order of magnitude as 
the mean distance between objects of a new phase r mean . Really, when 5r 
attains imaginary value r mean it means that the process of nucleation stops. 
So, namely the values at this moment are the final ones and the main ones. 
Then 

Tmean ~ "J" 

The same will be valid for all other first order transitions. 

Consider the opposite transition from the liquid to the vapor phase. We 
shall see the opposite inequality. But we know that compressibility of liquid 
is very law. So, the relative super streching is rather low. It means that the 
vapor phase in a final state will occupy some rather small volume. And again 
we see that the mean distance between the objects of a new phase is many 
times greater than the size of objects. Again we get the same result. 

We see that the following fact takes place: The relative final volume of a 
new phase will be small. 

To see that (^) takes place we needn't to consider the real final volume 
of a new phase after phase transition but can take the values at the end of 
nucleation. Because the relative quantity of surplus substance condensed in a 
new phase has an order T _1 (at least under the collective regime of substance 
consumption 0) we see that the relative volume occupied by a new phase is 
limited from above by a value of an order So 

r ~ i?r 1/3 

' mean 1 LJ - 

and the required property is established. 

So, we see the required property for moderate effects of depleted zones. 
When the effect of depletion zones is strong we can use at r ~ (2 4- 3)R 
stationary solution and see that 5r > (5-i-7)R. So, here the required property 
is also observed which completes justification. 

Then we can use profiles obtained on the formalism of Green function. 
We have § 

»(oo) - n(r) = J q ^ _ ^(-^—^ (6) 

A = (4vr) 3 / 2 (^) 1 / 2 (Cn 0O D) 3 / 2 (7) 
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and C is a supersaturation. 

This result was obtained for transition from vapor phase into a liquid 
phase but there is no difference because the kinetic mechanismremains the 
same. 

One can see that expressions (D, © and (|), © are absolutely different 
in analytical structure. The same is valid when we compare ([I]), (0) and (^|), 
(|7|) numerically. 

On the base of profile we can easily get all results of the nucleation process. 
The rate of nucleation can be found according to (§). Then we have 

J(r) = I(r = oo) exp(rn(0) — n(oo)— er/c(— ^=)/n(oo)) (8) 

for analytical solution with fixed boundary of embryo and 

7(r) = J(r = oo)exp(-r^ %{D ^ t ^ x)f/2 ex P (-^-^)dx/n(oo)) 

(9) 

for solution on the base of Green functions. Note that the analogous trans- 
formation in || was done with partial shift of coordinate r — R to r which is 
wrong and beside this isn't necessary. 

We see that these approaches give absolutely different results. 

The next step is an obvious remark that the probability dp that an embryo 
appears during elementary interval dt at the distance from r up to r + dr 
from the center of the already appeared embryo is 

dp = dtA7rr 2 I(r, t) 

Then the probability dP that the embryo appears in the layer from R up to 
r is 

dP = dt 47Tr' 2 I(r',t)dr' (10) 
Jr 

This expression differs from Q where the following expression was pre- 
sented 



dP = dt f 4nr 2 I{r',t)dr' 

Jr 



The last expression is wrong and we shall follow (|T0|). 

The next step used in [Q is to come to an integral value P which is the 
probability that in a sphere of radius r around already existing embryo there 
will be no appearance of a new embryo. Certainly we can start from (|10"D 
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and integrate it. But it is necessarfy to take into account that the zone of 
depletion will grow according to the absolutely precise law. This law was 
established in ||. But here we supposeQthat one can act in another manner 
and write 

P = exp(- [ [ r 4irr' 2 I(r',t')dr'dt') (11) 

JO JR 

The last expression analogous to those which forms the base for further 
constructions inQ || is rather doubtful. The problem appears because now 
we are coming to characteristic of a whole process of nucleation (earlier there 
was a profile around a solitary object of a new phase). 

For a "solitary subsystem" the derivation of (|TT|) is evident: Let 

dp = dtl 

be a probability that during an elementary interval dt there will be appear- 
ance of a new embryo. The value of dt is rather small. Then dp is proportional 
to dt and to the rate of nucleation / (here the sense of the nucletaion rate as 
a probability is used). Then the probability for the absence of appearance is 

dp' = 1 — dp = 1 — dtl 

If we have two elementary intervals dt\ and dti which don't overlap then the 
probability of the absence of embryo appearance is 

dp'\dp' 2 = (1 — I\dti){l — Iidt<i) 

The same is valid for arbitrary number of intervals. 

To fulfill multiplication in the r.h.s. one can present dp 1 as 

dp' = (1 — Idt) = exp(-Idt) 

Then the total probabiblity of the absence of appearance is 

P = IT d ti = ex P(- E dt ili) = exp(- J I{t)dt) (12) 

i i 

and we come to (O). 
lr To be close to Q. 

2 In B an intergal over r is from up to r. 
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The real problem is how to take into account the overlapping of density 
profiles initiated by different embryos. This problem was solved in 0. The 
exponential form ( [TT]) for P is some rather arbitrary interpretation of the 
value "the free volume" used in ||. Now we shall discuss this interpretation 
in frames of the approximation of solitary droplet. 

The overlapping of exhausted regions (ER) requires to use instead of (|8]) - 
(||)another approximation where instead of one profile there is a superposition 
of profiles. Certainly it is impossible to calculate this superposition precisely. 
Then one will come to some models analogous to formulated in ||. 

One can directly use the results from |2[] here. When the free volume |2j 
is calculated as function of time one can determine the total nucleation rate 
as a ratio of the free volume to the whole volume of a system and use then 

P). 

The real problem is what consequences can we make from a knowledge 
of P. On one hand it seems that the knowledge of P solves al problems in 
knetics of the global nucleation stage. Really, having presented P in a form 

P = exp(-L(r, t)) 

where L(r,t) is some expression one can approximately estimate the average 
radius of a sphere where there will be no appearance by relation 

L(r,t -> oo)| r=ro = 1 

This construction analogous to |4| needs two remarks. 

This estimate makes no difference between embryos appeared earlier or 
later during the nucleation stage. It means that this approximation supposes 
that the spectrum of droplets sizes is monodisperce one. It isn't correct 
assumption and can be treated only as a rough estimate. 

It is known that sometimes the overlapping of profiles is important and 
makes the main contribution in nucleation kinetics. This phenomena can 
occur when the substance exchange regime is goimg to be a free molecular 
one and when the long tails of profiles are important. The first situation can 
not be realized here. The last situation was a metter of separate consideration 
in@. 

In [0 an integral definition of a boundary of ER allowed to take into 
account the situation of long tails. Here this situation can not be considered, 
because the level type definition is used. It means that the value of r is 
determines as a value when L attains some level (= 1). The definition of 
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level type can not take into account long tails and it is more convenient 
to use approach from 0. Here this approach is considered because it is 
analogous to ||]. 

When the value of tq is determined it is easy to estimate the total average 
number of droplets N as 

3 

47r' 

(the coefficient can be omitted, this depends on interpretation of r ). 

To fulfill concrete calculations one can use the following approach. Note 
that P can be in any case presented as 



N = — r^ 3 



poo pr 

exp(— / dt' / dr'Aur' 2 exp(/(r' ,t')) 

JO Jri 



The lower limit 77 of integration has to be put 77 = R for models with 
explicit boundary of an embryo and for models with Green function 77 = 0. 
But because 5r 3> R the effect of the lower limit will be small and we can 
put 77 = in 11 situations. Certainly, for P due to the exponential form 
the effect is essential, but when we are calculating the mean distances the 
essential dependence on 77 disappears due to 5r ^> R. In the last equation / 
is some function with explicit form given by solution of diffusion equation. 

One can see that in f(r,t) two variables has to appear in combination 
(3 = t/r 2 or f3' = t/{r — R) 2 . The dependence on t and r via f3 is the main one 
(in special regimes of mother phase consumption there may be dependence on 
t, r in another combination, but this dependence will be much more smooth 
than via 0). Because SR R the difference in substitution (3' instead of j3 
will be small (not in P, but on mean distances and times). 

Then in integration f dt' we have to substitute dt' by d(3' which gives 

r r 

P = exp(— / dr'Anr' 4 const 
Jo 

where const comes from 

const = I d(3'p' 2 exp(f(r',t')) 

because 

It leads to 







OO POO 

dt' = r' 2 / dp' 
Jo 



P = exp(-r 5 7 ) (13) 
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where 7 = Airconst/5 is some constant. 
One can show that _ 

poo 

/ dxexp(— x 5 ) ~ 1 
Jo 

This approximate equality is very important. Namely this equality allows to 
determine the mean distance between droplets as iV~ 1//3 , or as r or according 
to the sense of P. Three ways are available. They have to give similar results. 
This fact is ensured by this approximate coincidence. Then the third way 
gives 

roo (J^ roo 

f = J Q drPr = ^ JE J o daexp(-a 5 / 2 ) 

where Co is the normalizing factor of distribution P. One can calculate Co 
according to 

/•OO 

Co / <irexp(— 7r 5 ) = 1 
Jo 



and 



^1/5 

Co 



Jq 00 <ir exp(— x 5 ) 
Therfl 

^ r> _ 1 J °°daexp{-a 5 / 2 ) 



2 7 V5 / oo da;exp(-a; 5 ) 
Because 

/ daexp(-a 5/2 ^ 
Jo 



POO 

\ dxexp(— x 5 ) 
Jo 



we see that f is two times smaller than the mean distance between embryos. 
It is clear because the mean distance between embryos is two mean distances 
until the boundary of exhaustion zone r . So, we observe a coincidence. 
Certainly, one can calculate integrals more precisely. 

The average time of waiting fro the appearance of new embryo near the 
already existing one is directly and elementary connected with a mean dis- 
tance between embryos and we needn't to calculate it separately^. 

3 This value is calculated in Q in a wrong way. 

4 Tlie integral in expression for this value is calculated in H in a wrong way 
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As for calculation a quadratic mean deviation from the average distance 
it is givenQ by 

< (r- < r >) 2 >=< r > 2 - 1) 

where _ 

Ii = dxx 1 exp(— a; 5 ) 
Jo 

But we have to stress that this result has practically no meaning because 
as it is stated in |7j the effect of interaction leads to decrease of fluctuation 
of the total number of droplets. The right numerical value of this effect is 
given in ||. Here we have absolutely no interaction between droplets and 
this result has to be seriously changed due to interaction. An example of 
such account is given in ||. 

To close this question we shall construct now more precise solutions of 
diffusion equation. They have to be used in (|J) and then in ( [111 ) 

We shall start with the following model. There exists a more general 
solution p. 

For a domain R < r < oo where 

n = f(r) at t = 

and 

n = g(t) at r = R 
the solution of diffusion equation 

dt dr 2 r dr 

(this the diffusion equation with spherical simmetry) is 

U= ^(0[exp(-i^)-ex P ( — 

(14) 

2R (r-R) 2 



where 



R 



2VDt 



5 Result in B is wrong. 
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This solution allows to formulate the following approximation. Let R be 
constant 

R = R(t) 

at the current moment when we want to know the density profile. But 
the boundary condition n b at R(t) in the previous moment t — t' will be 
recalculated according to the stationary solution: 

, . n(oo) — n(0) s . 
n » = n (°°) - R(t-t>) R(t) (15) 

Really the rate of embryo growth is so small that the stationary distribu- 
tion can be regarded as practically precise one. So, this approximation is 
practically precise and takes into acound the variation of mother phase con- 
sumption by an embryo. 

Note that diffusion equation 

dw £)\® 2w 2<9u> 
dt dr 2 r dr 

in the system with spherical symmetry can be transformed by substitution 

u(r, t) = rw(r, t) 

to an ordinary one dimensional diffusion equation 

du d 2 u 
dt dr 2 

We shall consider this equation in future. 

Then we can make a shift x = r — R(0). The problem 

w = f(x) 

at t = 

w = g(t) 

at x = has solution 

w(x,t) = —^= /°°[exp(- (a: ~^ 2 ) - exp(- (z "t/ )2 )]/(0<% 
l^f^Wth 1 V ADt ' FK ADt nJy ^' s 



(16) 



x f* f x 2 g{j) 
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It is also possible to go to variable 

p = r- R(t) 

where R(t) = At 1 / 2 is a boundary of an embryo. Then in the diffusion 
equation beside 

~ 1 du d 2 u 
QU = DM ~ W 

there appear some additional terms. Having considered Q as a main operator 
we come to an iteration procedure where all other terms are assumed to be 
known in a previous approximation (in zero approximation they are omitted). 
So, at every step of iterations we have to solve a problem 

w = f(x) 

at t = 

w = g(t) 

at p = R(0) for equation 

1 du d 2 u 

may be with renormalized D. Here $' is a known (at the previous step) 
function. After the mentioned shift x = p — R(0) we came to a problem 

w = f(x) 

at t = 

w = g(t) 

at p = R{0) for equation 

1 du d 2 u 

with known $. The solution of the last problem is 

u= / G(x,U)f(0^ + 
Jo 

x [* r x 2 g(r) 



2^V 4 J D(t-r) J (t-r) 3 / 2 

rt roc 

+ / / G(x,Z,t-T)$(Z,T)d£dT 

Jo Jo 
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where 



1 



[exp( 



(x-0 2 



) - exp( 



(x + 2 



)] 



ADt 



ADt 



The last relation solves the problem to construct iterations on every step. 
Thus, we come to solution of diffusion equation which has to be used in our 
previous constructions. 

1 Numerical results 

We have calculated the worst situation when v v /v p = 10 for vapor-liquid 
transition. Here we have taken G — 13 which is the minimal value for 
macroscopic description of the critical embryo, when v v /v\ and G increase 
the accuracy of approach (H) and partially the accuracy of approach fll5| . 
The accuracy of ([]]) decreases because the error in (|l|) depends mainly on the 
property of relative stationarity which has no connection with values of Y. 
It isn't so evident that the property v v /vi 3> 1 also has no influence on the 
error of ([!]) because here only the relative variation of radius is important. 
But when V increases it means that the main important region of solution 
now corresponds to the asymptote where the quakitative behavior of (|l]) is 
wrong. So, we came to conclusion that it isn't possible to use ([!]). The 
error produced by ([]]) is seen already at T — 13. Figure 1 illustrates the 
behavior of different profiles in this situation. One can see four lines there. 
One line is apart from all other and it has the different asymptotic behavior. 
This is solution ([!]). Other lines are solutions of diffusion problem and some 
approximations. One of them (the highest) is thick because actually there 
are two formally coinciding solutions: precise solution and approximation 
flT5|). The coincidence is formal. In the middle of this group one can find 
numerical solution with prescribed law of droplet radius growth R = At 1 / 2 
(it is obtained by linear extrapolation of results of precise solution). The 
difference between these curve and precise self consistent solution is caused 
by numerical errors and corrections on radius growth, the lowest curve in this 
group is approximation (|6|). One can see that both ([]) and ([X5|) give good 
results, for natural situations with big values of Y and v v jv\ the accuracy 
will even increase. 
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Figure 1. 

Different approximations for the form of density profile. 
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